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We investigate the finite-temperature properties of attractive three-component (col- 
ors) fermionic atoms in optical lattices using a self-energy functional approach. As the 
strength of the attractive interaction increases in the low temperature region, a second- 
order transition occurs from a Fermi liquid to a color superfiuid (CSF). In the strong 
attractive region, a first-order transition occurs from a CSF to a trionic state. In the high 
temperature region, a crossover between a Fermi liquid and a trionic state is observed 
with increasing the strength of the attractive interaction. The crossover region for fixed 
temperature is almost independent of filling. 
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1. Introduction 

Recent significant progress in cold fermionic atoms loaded in optical lattices pro- 
vides us with a way to investigate the fascinating aspects of quantum many-body 
effects. Novel phenomena have been studied both experimentally and theoretically. 
A balanced mixture of three dif ferent hyperfine states of attractive 6 Li fermionic 
atoms was succeeded in creating The observed three-body loss was discussed 
in terms of Efimov states B B> | 6 | 7 | 8 | Q n ^ e Q-j^er hand, theoretical studies for attrac- 
tive three-component (col or degrees of fr eedom) fermionic atoms in optical lattices 
were performed intensively ^^^^^ l It was shown that a color superfiuid (CSF) 
emerges for the small attractive interaction, where atoms with two of three colors 
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form the Cooper pairs and the third ones remain a Fermi liquid. As the attrac- 
tive interaction bec omes stro ng, a quantum phase transition from the CSF to the 
trionic state occurs H HUSS j n th e trionic state, singlet bound states of three dif- 
ferent color atoms are formed. Despite the high tunability of the optical lattice 
systems, it is difficult to cool the fermionic system because of the Pauli principle. 
It is thus desirable to investigate finite-temperature properties of the attractive 
three-component fermionic atoms in optical lattices. 

In a recent study, we showed that a crossover occurs between the Fermi liquid 
and the trionic state at finite temperatures Furthermore, the crossover region 
in the phase diagram concerning temperature and interaction hardly depends on 
temperature at half filling. In real systems, a harmonic confinement potential exists, 
which induces a spatial change in filling. To observe the trionic state in experiments, 
the effects of temperature and filling change are indispensably investigated. 

In this paper, we investigate the CSF and trionic state of attractive three- 
component fermionic atoms in optical lattices at finite temperatures. We discuss 
the stability of the trionic state against filling change at finite temperatures. 

2. Model and Method 

The system is assumed to keep a balanced population of fermionic atoms with three 
colors. The fermionic optical lattice systems can be well described by the following 
Hubbard-type Hamiltonian 031 

n = -t 4a c ja n ia ri0 ~^2fi a n ia , (1) 

(i,j),a i,a^P i,a 

where Ci a annihilates a fermionic atom with color a{— 1, 2, 3) at the ith lattice site, 
ni a = cl a Ci a , the subscript indicates the sum over the nearest-neighbor sites, 
and — {/(< 0) is the on-site attractive interaction between two atoms with different 
colors. 

Using a self-energy functional approach (SFA) I 15 | 16 |^ we elucidate character- 
istics of the CSF, trionic state, and Fermi liquid, and study the phase transition 
and crossover between them. The SFA allows us to deal efficiently with zero- and 
finite-t emper ature properties concerning the phase transition driven by correlation 
effects indeed, precise results related to the Mott transition in correlated 

electron systems were obtained as regards thermodynamic quantities, excitation 
spectra, and phase diagrams. This method is based on the Luttinger-Ward varia- 
tional method. The framework of this variational approach enables us to introduce 
a proper reference system which has to include the same interacting Hamiltonian 
as that of the original Hamiltonian. Here, we apply the two-site Anderson impurity 
model as the reference system. Para meters of the reference system are determined 
by the variational method. In the sfa™, the grand potential VL can be expressed 
as a function of the self-energy of the reference system Et ■ From the condition that 
the derivative of the grand potential with respect to the variational parameters be- 
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comes zero, dVl/dt — (dUt / dt){d£l / dUt) — 0, we obtain the reference self-energy 
St which properly describes the original correlated system, where t is the parameter 
matrix of the nonintcracting part of the reference Hamiltonian. 

We use a semicircular density of states obtained for the infinite-dimensional 
Bethe lattice. We calculate the CSF order parameter $ = (cjjcjj), the quasiparticle 
weight Z = [1 — d'St(oj)/duj}\~L , which is independent of color, and the entropy 
per site S/L = —d(Q/L)/dT, where L is the number of lattice sites. We further 
calculate the expectation value of the trion number per site N t — (N) and the 
local susceptibility of trion-number fluctuations defined by \ = J e lUT (: N(t) :: 
Ni(Q) :}dr\ u=a , where : N(t) := N(t) - N t with N = rmn^n^. Note that $, N t , 
and x are independent of i and color in the SFA procedure. We do not consider 
the possibility of the trionic density- wave state induced by spatial fluctuations. The 
hopping integral t is used in units of energy. 

3. Results and Discussions 




Fig. 1. \U\/t dependence of (a) the CSF order parameter <!> and (b) the quasiparticle weight Z at 
T = 0. 

We begin our discussions with zero-temperature properties. The chemical po- 
tential is set at fi = U so that particle-hole symmetry can be satisfied and that 
a balanced population of each color atoms can be achieved. Figure [T] shows the 
CSF order parameter <I> and the quasiparticle weight Z as a function of \U\/t. As 
\U\/t is increased, $ first increases and vanishes discontinuously at U/t — —1.94. 
The quasiparticle weight Z decreases significantly and then drops discontinuously 
to zero at the same U/t — —1.94. The results indicate that the CSF disappears as a 
result of a first-order Mott transition driven by the renormalization effects. In Fig. 
HJa), we show the expectation value of the trion number Nt as a function of \U\/t. 
At half filling, the maximum of Nt is 0.5. As \U\/t increases, Nt increases and jumps 
to ~ 0.5 at the same critical U/t. We thus conclude that the first-order CSF-trion 
quantum phase transition occurs at U/t = —1.94. These results are in contrast with 
those in one-dimensional systems, where the CSF disappears at zero temperature 
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Fig. 2. \U\/t dependence of (a) the expectation value of the trion number Nt and (b) the local 
susceptibility of trion-number fluctuations \ at T = 0. 

because of strong quantum fluctuatio ns and the t hree-particle density-wave states 
caused by spatial fluctuations appear | 17 | 18 | 19 |20| 




Fig. 3. Grand potentials for (a) U/t = —1.3 and (b) U/t = —1.9 as a function of temperature. FL 
means Fermi liquid. 



We investigate the phase transition and crossover at finite temperatures. In Fig. 
G2a) and (b), we show the grand potentials for U/t = —1.3 and U/t = —1.9 as a 
function of temperature. In the low temperature region for U/t = —1.3, the grand 
potential of the CSF is lower than that of the Fermi liquid. When temperature 
is increased, both grand potentials agree continuously. Therefore, for U/t = — 1.3 
the second-order phase transition occurs from the CSF to the Fermi liquid with 
increasing temperature. For U/t = —1.9, the grand potential of the CSF is lower 
than that of the trionic state in the low temperature region. When temperature is 
increased, both grand potentials cross, indicating that a first-order phase transition 
occurs from the CSF to the trionic state. 

In Fig. Eta), the contour lines of the entropy are shown in the T — \U\ plane. 
For \U\/t < 1.7, we find the cusp in each line. This is a manifestation of the second- 
order Fermi liquid-CSF transition. For T/t > 0.08 around \U\/t = 1.7, we find the 
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Fig. 4. (a) Contour lines of the entropy S/L in the T-|[/| plane, (b) Phase diagram of the isotropic 
attractive three-component Hubbard model at half filling. FL means Fermi liquid. 



large entropy region where S/L is above 0.8. When \U\/t is increased, the entropy 
decreases continuously to S/L ~ In 2 ~ 0.693. Accordingly, a crossover occurs 
from the Fermi liquid to the trionic state. We summarize the results calculated for 
various parameters in the phase diagram for the Fermi liquid, CSF, and trionic 
state shown in Fig. lUJb). The second-order phase transition occurs between the 
Fermi liquid and the CSF, and the first-order phase transition occurs between the 
CSF and the trionic state. There occurs a crossover between the Fermi liquid and 
the trionic state. The crossover region seems to be \U\/t ~ 2. This is because the 
crossover occurs when the trionic excitation gap (2|/7|) and the bandwidth (At) 
become comparable to each other. Because of thermal fluctuations, the crossover 
region spreads gradually with increasing T/t. In the phase diagram, the amount 
of the CSF order parameter $ is also shown. As \U\/t is increased, $ increases 
gradually from zero and vanishes suddenly from the value close to the maximum. 
This result suggests that the CSF is suppressed by the trion formation. In a recent 
paper, we demonstrated that, when the anisotropy of the attractive interactions are 
included, the CSF phase extends actually towards the large \U\/t and T/t region 

[To] 

because of the suppression of the trion formation 

We now discuss the origin of S/L = In 2 in the trionic state at half filling. In 
Fig. [2jb), we show the local susceptibility of trion- number fluctuations % at T = 0. 
When \U\/t approaches 1.94 from below, \ shows divergence. This divergence is 
caused by the formation of a localized trion. Within the SEA, the trionic state can 
be regarded as the Mott insulating state as shown in Fig. [Tib), where the localized 
trion is formed at a site. This localized trion induces two degrees of freedom whether 
or not the trion is found at a site. Therefore, the entropy in the trionic state is In 2 
in T = at half filling. A similar behavior is seen in the spin susceptibility of 
strongly correlated electron systems. When the two-component repulsive Hubbard 
model is treated by a dynamical mean-field theory, the local spin susceptibility 
diverges with the system approaching the Mott transition point from the weakly 
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interacting regional This is because a localized 5 = 1/2 free spin is formed in 
the Mott insulating state in the infinite dimensional system, yielding the residual 
entropy In 2. 

The residual entropy is caused by neglecting spatial fluctuations in the present 
SFA calculation. The present calculation is valid for the infinite dimensional system 
or for the high temperature region, where spatial fluctuations are suppressed. By 
contrast, in the finite dimensional syst em, spatial fluctuations cause sometimes 
relevant effects at low temperatures L^M As mentioned before, the trionic density- 
wave state is expected to appear in low temperatures at half filling. Furthermore, 
the trion hopping is expected to be induced by spatial fluctuations. The effective 
hopping integral and the effective nearest-neighbor repu lsio n can be derived as t c g ~ 
t 3 /U 2 and Uos — t 2 /\U\ for large \U\/t, respectively^. Therefore, the relation 
t c ff < U e g is satisfied, suggesting strongly correlated trions. It is thus expected that 
in the finite dimensional systems the trionic density-wave state appears close to 
half filling in low temperatures and the strongly renormalized trionic Fermi liquid 
appears in T < t e s away from half filling. For low filling, the effective repulsion 
may become irrelevant and the strong renormalization effects may be suppressed. 
The investigation of the present system including the effects of spatial fluctuations 
is our future study. 





Fig. 5. Entropy S/L at T/t = 0.08 as a function of \U\/t (a) for filling < n < 1, where n = 0.5 
denotes half filling, and (b) for fixed Sfi. 



Finally, we investigate the crossover between the Fermi liquid and the trionic 
state for various filling with a balanced population of each color atoms. We calculate 
the entropy for < n < 1, where n is filling and n = 0.5 denotes half filling. To 
obtain the accurate results, we set at T/t = 0.08 and neglect the possibility of the 
CSF. The results are shown in Fig. [5] as a function of \U\/t. In Fig. [5ja), we find 
that the maximum of the entropy S/L appears in 0.1 < n < 0.9. When \U\/t is 
increased in this region, S/L for fixed n increases gradually and takes a maximum 
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Table 1. \U\/t and n of the maximum shown in Fig.(5fb). 





\u\/t 


n 





1.7 


0.5 


0.1 


1.7 


0.64 


0.2 


1.7 


0.77 


0.3 


l.G 


0.86 



in 1.6 < \U\/t < 1.7. After the maximum, S/L decreases gradually, which indicates 
a crossover between the Fermi liquid and the trionic state. The results demonstrate 
that the crossover region appears in 1.6 < \U\/t < 1.7 irrespective of filling for 
fixed temperature. In Fig. EJb), we show S/L for fixed it, where 5fi = fi — U is 
the shift of the chemical potential from that for half filling. For half filling = 
0), S/L approaches In 2 with increasing \U\/t. For 5/i ^ 0, S/L decreases after a 
maximum. This is because particle-hole symmetry of the trion is broken and filling 
decreases with increasing \U\/t for fixed 8[i ^ 0. In Tab. 1, we show \U\/t and n 
of the maximum shown in Fig. [5th). We find that \U\/t of the maximum is almost 
independent of filling. In other words, the attractive interaction of the crossover for 
fixed temperature is almost independent of filling. 

In real systems, a harmonic confinement potential exists, which yields a spatial 
change in filling. The center of the confinement potential is so dense. The results 
in Fig. [SJa) effectively presents the effects of a confinement potential. Decrease in 
n denotes that the observed position becomes distant from the center of the con- 
finement potential. It is expected that the trionic state appears in the rather wide 
spatial range corresponding to 0.1 < n < 0.9. We have shown that the crossover 
between the Fermi liquid and the trionic state is almost independent of filling. We 
expect that in real experiments the crossover occurs at nearly the same interac- 
tion irrespective of spatial change in filling. In addition, as shown in Fig.[4jb), the 
crossover occurs at \U\/t ~ 1.5 — 2.0 irrespective of temperature. These findings are 
useful for detecting the trionic state and the crossover. 

The large entropy ~ In 2 of the trionic state close to half filling suggests the 
possibility of the ad iabatic cooling during the lattice loading caused by the Pomer- 
anchuk effect BUM] However, the entropy of the trionic state decreases away from 
half filling. It is interesting to investigate the competition between the adiabatic 
cooling and the reduction of the entropy caused by the confinement potential. This 
issue constitutes our future study. 

4. Summary 

We have investigated the attractive three-component fermionic atoms in optical 
lattices at zero and finite temperatures. We have shown that the attractive interac- 
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tion of the crossover region between the Fermi liquid and the trionic state for given 
temperature is almost independent of filling. The present results are complementary 
to our previous results ^1 We hope that our results will contribute to the study of 
the exotic phases of multicomponent cold fermionic atoms in optical lattices. 
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